home *** CD-ROM | disk | FTP | other *** search
/ Windows News 2010 Summer - Disc 1 / WN_Ete2010_CD1.iso / Onglet5 / Weezo / Weezo setup.exe / {code_appDir} / www / js / BigInt-src.js < prev    next >
Text File  |  2010-05-19  |  16KB  |  630 lines
  1. // BigInt, a suite of routines for performing multiple-precision arithmetic in
  2. // JavaScript.
  3. //
  4. // Copyright 1998-2005 David Shapiro.
  5. //
  6. // You may use, re-use, abuse,
  7. // copy, and modify this code to your liking, but please keep this header.
  8. // Thanks!
  9. //
  10. // Dave Shapiro
  11. // dave@ohdave.com
  12.  
  13. // IMPORTANT THING: Be sure to set maxDigits according to your precision
  14. // needs. Use the setMaxDigits() function to do this. See comments below.
  15. //
  16. // Tweaked by Ian Bunning
  17. // Alterations:
  18. // Fix bug in function biFromHex(s) to allow
  19. // parsing of strings of length != 0 (mod 4)
  20.  
  21. // Changes made by Dave Shapiro as of 12/30/2004:
  22. //
  23. // The BigInt() constructor doesn't take a string anymore. If you want to
  24. // create a BigInt from a string, use biFromDecimal() for base-10
  25. // representations, biFromHex() for base-16 representations, or
  26. // biFromString() for base-2-to-36 representations.
  27. //
  28. // biFromArray() has been removed. Use biCopy() instead, passing a BigInt
  29. // instead of an array.
  30. //
  31. // The BigInt() constructor now only constructs a zeroed-out array.
  32. // Alternatively, if you pass <true>, it won't construct any array. See the
  33. // biCopy() method for an example of this.
  34. //
  35. // Be sure to set maxDigits depending on your precision needs. The default
  36. // zeroed-out array ZERO_ARRAY is constructed inside the setMaxDigits()
  37. // function. So use this function to set the variable. DON'T JUST SET THE
  38. // VALUE. USE THE FUNCTION.
  39. //
  40. // ZERO_ARRAY exists to hopefully speed up construction of BigInts(). By
  41. // precalculating the zero array, we can just use slice(0) to make copies of
  42. // it. Presumably this calls faster native code, as opposed to setting the
  43. // elements one at a time. I have not done any timing tests to verify this
  44. // claim.
  45.  
  46. // Max number = 10^16 - 2 = 9999999999999998;
  47. //               2^53     = 9007199254740992;
  48.  
  49. var biRadixBase = 2;
  50. var biRadixBits = 16;
  51. var bitsPerDigit = biRadixBits;
  52. var biRadix = 1 << 16; // = 2^16 = 65536
  53. var biHalfRadix = biRadix >>> 1;
  54. var biRadixSquared = biRadix * biRadix;
  55. var maxDigitVal = biRadix - 1;
  56. var maxInteger = 9999999999999998; 
  57.  
  58. // maxDigits:
  59. // Change this to accommodate your largest number size. Use setMaxDigits()
  60. // to change it!
  61. //
  62. // In general, if you're working with numbers of size N bits, you'll need 2*N
  63. // bits of storage. Each digit holds 16 bits. So, a 1024-bit key will need
  64. //
  65. // 1024 * 2 / 16 = 128 digits of storage.
  66. //
  67.  
  68. var maxDigits;
  69. var ZERO_ARRAY;
  70. var bigZero, bigOne;
  71.  
  72. function setMaxDigits(value)
  73. {
  74.     maxDigits = value;
  75.     ZERO_ARRAY = new Array(maxDigits);
  76.     for (var iza = 0; iza < ZERO_ARRAY.length; iza++) ZERO_ARRAY[iza] = 0;
  77.     bigZero = new BigInt();
  78.     bigOne = new BigInt();
  79.     bigOne.digits[0] = 1;
  80. }
  81.  
  82. setMaxDigits(20);
  83.  
  84. // The maximum number of digits in base 10 you can convert to an
  85. // integer without JavaScript throwing up on you.
  86. var dpl10 = 15;
  87. // lr10 = 10 ^ dpl10
  88. var lr10 = biFromNumber(1000000000000000);
  89.  
  90. function BigInt(flag)
  91. {
  92.     if (typeof flag == "boolean" && flag == true) {
  93.         this.digits = null;
  94.     }
  95.     else {
  96.         this.digits = ZERO_ARRAY.slice(0);
  97.     }
  98.     this.isNeg = false;
  99. }
  100.  
  101. function biFromDecimal(s)
  102. {
  103.     var isNeg = s.charAt(0) == '-';
  104.     var i = isNeg ? 1 : 0;
  105.     var result;
  106.     // Skip leading zeros.
  107.     while (i < s.length && s.charAt(i) == '0') ++i;
  108.     if (i == s.length) {
  109.         result = new BigInt();
  110.     }
  111.     else {
  112.         var digitCount = s.length - i;
  113.         var fgl = digitCount % dpl10;
  114.         if (fgl == 0) fgl = dpl10;
  115.         result = biFromNumber(Number(s.substr(i, fgl)));
  116.         i += fgl;
  117.         while (i < s.length) {
  118.             result = biAdd(biMultiply(result, lr10),
  119.                            biFromNumber(Number(s.substr(i, dpl10))));
  120.             i += dpl10;
  121.         }
  122.         result.isNeg = isNeg;
  123.     }
  124.     return result;
  125. }
  126.  
  127. function biCopy(bi)
  128. {
  129.     var result = new BigInt(true);
  130.     result.digits = bi.digits.slice(0);
  131.     result.isNeg = bi.isNeg;
  132.     return result;
  133. }
  134.  
  135. function biFromNumber(i)
  136. {
  137.     var result = new BigInt();
  138.     result.isNeg = i < 0;
  139.     i = Math.abs(i);
  140.     var j = 0;
  141.     while (i > 0) {
  142.         result.digits[j++] = i & maxDigitVal;
  143.         i = Math.floor(i / biRadix);
  144.     }
  145.     return result;
  146. }
  147.  
  148. function reverseStr(s)
  149. {
  150.     var result = "";
  151.     for (var i = s.length - 1; i > -1; --i) {
  152.         result += s.charAt(i);
  153.     }
  154.     return result;
  155. }
  156.  
  157. var hexatrigesimalToChar = new Array(
  158.  '0', '1', '2', '3', '4', '5', '6', '7', '8', '9',
  159.  'a', 'b', 'c', 'd', 'e', 'f', 'g', 'h', 'i', 'j',
  160.  'k', 'l', 'm', 'n', 'o', 'p', 'q', 'r', 's', 't',
  161.  'u', 'v', 'w', 'x', 'y', 'z'
  162. );
  163.  
  164. function biToString(x, radix)
  165.     // 2 <= radix <= 36
  166. {
  167.     var b = new BigInt();
  168.     b.digits[0] = radix;
  169.     var qr = biDivideModulo(x, b);
  170.     var result = hexatrigesimalToChar[qr[1].digits[0]];
  171.     while (biCompare(qr[0], bigZero) == 1) {
  172.         qr = biDivideModulo(qr[0], b);
  173.         digit = qr[1].digits[0];
  174.         result += hexatrigesimalToChar[qr[1].digits[0]];
  175.     }
  176.     return (x.isNeg ? "-" : "") + reverseStr(result);
  177. }
  178.  
  179. function biToDecimal(x)
  180. {
  181.     var b = new BigInt();
  182.     b.digits[0] = 10;
  183.     var qr = biDivideModulo(x, b);
  184.     var result = String(qr[1].digits[0]);
  185.     while (biCompare(qr[0], bigZero) == 1) {
  186.         qr = biDivideModulo(qr[0], b);
  187.         result += String(qr[1].digits[0]);
  188.     }
  189.     return (x.isNeg ? "-" : "") + reverseStr(result);
  190. }
  191.  
  192. var hexToChar = new Array('0', '1', '2', '3', '4', '5', '6', '7', '8', '9',
  193.                           'a', 'b', 'c', 'd', 'e', 'f');
  194.  
  195. function digitToHex(n)
  196. {
  197.     var mask = 0xf;
  198.     var result = "";
  199.     for (i = 0; i < 4; ++i) {
  200.         result += hexToChar[n & mask];
  201.         n >>>= 4;
  202.     }
  203.     return reverseStr(result);
  204. }
  205.  
  206. function biToHex(x)
  207. {
  208.     var result = "";
  209.     var n = biHighIndex(x);
  210.     for (var i = biHighIndex(x); i > -1; --i) {
  211.         result += digitToHex(x.digits[i]);
  212.     }
  213.     return result;
  214. }
  215.  
  216. function charToHex(c)
  217. {
  218.     var ZERO = 48;
  219.     var NINE = ZERO + 9;
  220.     var littleA = 97;
  221.     var littleZ = littleA + 25;
  222.     var bigA = 65;
  223.     var bigZ = 65 + 25;
  224.     var result;
  225.  
  226.     if (c >= ZERO && c <= NINE) {
  227.         result = c - ZERO;
  228.     } else if (c >= bigA && c <= bigZ) {
  229.         result = 10 + c - bigA;
  230.     } else if (c >= littleA && c <= littleZ) {
  231.         result = 10 + c - littleA;
  232.     } else {
  233.         result = 0;
  234.     }
  235.     return result;
  236. }
  237.  
  238. function hexToDigit(s)
  239. {
  240.     var result = 0;
  241.     var sl = Math.min(s.length, 4);
  242.     for (var i = 0; i < sl; ++i) {
  243.         result <<= 4;
  244.         result |= charToHex(s.charCodeAt(i))
  245.     }
  246.     return result;
  247. }
  248.  
  249. function biFromHex(s)
  250. {
  251.     var result = new BigInt();
  252.     var sl = s.length;
  253.     for (var i = sl, j = 0; i > 0; i -= 4, ++j) {
  254.         result.digits[j] = hexToDigit(s.substr(Math.max(i - 4, 0), Math.min(i, 4)));
  255.     }
  256.     return result;
  257. }
  258.  
  259. function biFromString(s, radix)
  260. {
  261.     var isNeg = s.charAt(0) == '-';
  262.     var istop = isNeg ? 1 : 0;
  263.     var result = new BigInt();
  264.     var place = new BigInt();
  265.     place.digits[0] = 1; // radix^0
  266.     for (var i = s.length - 1; i >= istop; i--) {
  267.         var c = s.charCodeAt(i);
  268.         var digit = charToHex(c);
  269.         var biDigit = biMultiplyDigit(place, digit);
  270.         result = biAdd(result, biDigit);
  271.         place = biMultiplyDigit(place, radix);
  272.     }
  273.     result.isNeg = isNeg;
  274.     return result;
  275. }
  276.  
  277. function biDump(b)
  278. {
  279.     return (b.isNeg ? "-" : "") + b.digits.join(" ");
  280. }
  281.  
  282. function biAdd(x, y)
  283. {
  284.     var result;
  285.  
  286.     if (x.isNeg != y.isNeg) {
  287.         y.isNeg = !y.isNeg;
  288.         result = biSubtract(x, y);
  289.         y.isNeg = !y.isNeg;
  290.     }
  291.     else {
  292.         result = new BigInt();
  293.         var c = 0;
  294.         var n;
  295.         for (var i = 0; i < x.digits.length; ++i) {
  296.             n = x.digits[i] + y.digits[i] + c;
  297.             result.digits[i] = n % biRadix;
  298.             c = Number(n >= biRadix);
  299.         }
  300.         result.isNeg = x.isNeg;
  301.     }
  302.     return result;
  303. }
  304.  
  305. function biSubtract(x, y)
  306. {
  307.     var result;
  308.     if (x.isNeg != y.isNeg) {
  309.         y.isNeg = !y.isNeg;
  310.         result = biAdd(x, y);
  311.         y.isNeg = !y.isNeg;
  312.     } else {
  313.         result = new BigInt();
  314.         var n, c;
  315.         c = 0;
  316.         for (var i = 0; i < x.digits.length; ++i) {
  317.             n = x.digits[i] - y.digits[i] + c;
  318.             result.digits[i] = n % biRadix;
  319.             // Stupid non-conforming modulus operation.
  320.             if (result.digits[i] < 0) result.digits[i] += biRadix;
  321.             c = 0 - Number(n < 0);
  322.         }
  323.         // Fix up the negative sign, if any.
  324.         if (c == -1) {
  325.             c = 0;
  326.             for (var i = 0; i < x.digits.length; ++i) {
  327.                 n = 0 - result.digits[i] + c;
  328.                 result.digits[i] = n % biRadix;
  329.                 // Stupid non-conforming modulus operation.
  330.                 if (result.digits[i] < 0) result.digits[i] += biRadix;
  331.                 c = 0 - Number(n < 0);
  332.             }
  333.             // Result is opposite sign of arguments.
  334.             result.isNeg = !x.isNeg;
  335.         } else {
  336.             // Result is same sign.
  337.             result.isNeg = x.isNeg;
  338.         }
  339.     }
  340.     return result;
  341. }
  342.  
  343. function biHighIndex(x)
  344. {
  345.     var result = x.digits.length - 1;
  346.     while (result > 0 && x.digits[result] == 0) --result;
  347.     return result;
  348. }
  349.  
  350. function biNumBits(x)
  351. {
  352.     var n = biHighIndex(x);
  353.     var d = x.digits[n];
  354.     var m = (n + 1) * bitsPerDigit;
  355.     var result;
  356.     for (result = m; result > m - bitsPerDigit; --result) {
  357.         if ((d & 0x8000) != 0) break;
  358.         d <<= 1;
  359.     }
  360.     return result;
  361. }
  362.  
  363. var inside_pow_mod = 0;
  364.  
  365. function biMultiply(x, y)
  366. {
  367.     var result = new BigInt();
  368.     var c;
  369.     var n = biHighIndex(x);
  370.     var t = biHighIndex(y);
  371.     var u, uv, k;
  372.  
  373.     for (var i = 0; i <= t; ++i) {
  374.         c = 0;
  375.         k = i;
  376.         for (j = 0; j <= n; ++j, ++k) {
  377.             uv = result.digits[k] + x.digits[j] * y.digits[i] + c;
  378.             result.digits[k] = uv & maxDigitVal;
  379.             c = uv >>> biRadixBits;
  380.             //c = Math.floor(uv / biRadix);
  381.         }
  382.         result.digits[i + n + 1] = c;
  383.     }
  384.     // Someone give me a logical xor, please.
  385.     result.isNeg = x.isNeg != y.isNeg;
  386.     return result;
  387. }
  388.  
  389. function biMultiplyDigit(x, y)
  390. {
  391.     var n, c, uv;
  392.  
  393.     result = new BigInt();
  394.     n = biHighIndex(x);
  395.     c = 0;
  396.     for (var j = 0; j <= n; ++j) {
  397.         uv = result.digits[j] + x.digits[j] * y + c;
  398.         result.digits[j] = uv & maxDigitVal;
  399.         c = uv >>> biRadixBits;
  400.         //c = Math.floor(uv / biRadix);
  401.     }
  402.     result.digits[1 + n] = c;
  403.     return result;
  404. }
  405.  
  406. function arrayCopy(src, srcStart, dest, destStart, n)
  407. {
  408.     var m = Math.min(srcStart + n, src.length);
  409.     for (var i = srcStart, j = destStart; i < m; ++i, ++j) {
  410.         dest[j] = src[i];
  411.     }
  412. }
  413.  
  414. var highBitMasks = new Array(0x0000, 0x8000, 0xC000, 0xE000, 0xF000, 0xF800,
  415.                              0xFC00, 0xFE00, 0xFF00, 0xFF80, 0xFFC0, 0xFFE0,
  416.                              0xFFF0, 0xFFF8, 0xFFFC, 0xFFFE, 0xFFFF);
  417.  
  418. function biShiftLeft(x, n)
  419. {
  420.     var digitCount = Math.floor(n / bitsPerDigit);
  421.     var result = new BigInt();
  422.     arrayCopy(x.digits, 0, result.digits, digitCount,
  423.               result.digits.length - digitCount);
  424.     var bits = n % bitsPerDigit;
  425.     var rightBits = bitsPerDigit - bits;
  426.     for (var i = result.digits.length - 1, i1 = i - 1; i > 0; --i, --i1) {
  427.         result.digits[i] = ((result.digits[i] << bits) & maxDigitVal) |
  428.                            ((result.digits[i1] & highBitMasks[bits]) >>>
  429.                             (rightBits));
  430.     }
  431.     result.digits[0] = ((result.digits[i] << bits) & maxDigitVal);
  432.     result.isNeg = x.isNeg;
  433.     return result;
  434. }
  435.  
  436. var lowBitMasks = new Array(0x0000, 0x0001, 0x0003, 0x0007, 0x000F, 0x001F,
  437.                             0x003F, 0x007F, 0x00FF, 0x01FF, 0x03FF, 0x07FF,
  438.                             0x0FFF, 0x1FFF, 0x3FFF, 0x7FFF, 0xFFFF);
  439.  
  440. function biShiftRight(x, n)
  441. {
  442.     var digitCount = Math.floor(n / bitsPerDigit);
  443.     var result = new BigInt();
  444.     arrayCopy(x.digits, digitCount, result.digits, 0,
  445.               x.digits.length - digitCount);
  446.     var bits = n % bitsPerDigit;
  447.     var leftBits = bitsPerDigit - bits;
  448.     for (var i = 0, i1 = i + 1; i < result.digits.length - 1; ++i, ++i1) {
  449.         result.digits[i] = (result.digits[i] >>> bits) |
  450.                            ((result.digits[i1] & lowBitMasks[bits]) << leftBits);
  451.     }
  452.     result.digits[result.digits.length - 1] >>>= bits;
  453.     result.isNeg = x.isNeg;
  454.     return result;
  455. }
  456.  
  457. function biMultiplyByRadixPower(x, n)
  458. {
  459.     var result = new BigInt();
  460.     arrayCopy(x.digits, 0, result.digits, n, result.digits.length - n);
  461.     return result;
  462. }
  463.  
  464. function biDivideByRadixPower(x, n)
  465. {
  466.     var result = new BigInt();
  467.     arrayCopy(x.digits, n, result.digits, 0, result.digits.length - n);
  468.     return result;
  469. }
  470.  
  471. function biModuloByRadixPower(x, n)
  472. {
  473.     var result = new BigInt();
  474.     arrayCopy(x.digits, 0, result.digits, 0, n);
  475.     return result;
  476. }
  477.  
  478. function biCompare(x, y)
  479. {
  480.     if (x.isNeg != y.isNeg) {
  481.         return 1 - 2 * Number(x.isNeg);
  482.     }
  483.     for (var i = x.digits.length - 1; i >= 0; --i) {
  484.         if (x.digits[i] != y.digits[i]) {
  485.             if (x.isNeg) {
  486.                 return 1 - 2 * Number(x.digits[i] > y.digits[i]);
  487.             } else {
  488.                 return 1 - 2 * Number(x.digits[i] < y.digits[i]);
  489.             }
  490.         }
  491.     }
  492.     return 0;
  493. }
  494.  
  495. function biDivideModulo(x, y)
  496. {
  497.     var nb = biNumBits(x);
  498.     var tb = biNumBits(y);
  499.     var origYIsNeg = y.isNeg;
  500.     var q, r;
  501.     if (nb < tb) {
  502.         // |x| < |y|
  503.         if (x.isNeg) {
  504.             q = biCopy(bigOne);
  505.             q.isNeg = !y.isNeg;
  506.             x.isNeg = false;
  507.             y.isNeg = false;
  508.             r = biSubtract(y, x);
  509.             // Restore signs, 'cause they're references.
  510.             x.isNeg = true;
  511.             y.isNeg = origYIsNeg;
  512.         } else {
  513.             q = new BigInt();
  514.             r = biCopy(x);
  515.         }
  516.         return new Array(q, r);
  517.     }
  518.  
  519.     q = new BigInt();
  520.     r = x;
  521.  
  522.     // Normalize Y.
  523.     var t = Math.ceil(tb / bitsPerDigit) - 1;
  524.     var lambda = 0;
  525.     while (y.digits[t] < biHalfRadix) {
  526.         y = biShiftLeft(y, 1);
  527.         ++lambda;
  528.         ++tb;
  529.         t = Math.ceil(tb / bitsPerDigit) - 1;
  530.     }
  531.     // Shift r over to keep the quotient constant. We'll shift the
  532.     // remainder back at the end.
  533.     r = biShiftLeft(r, lambda);
  534.     nb += lambda; // Update the bit count for x.
  535.     var n = Math.ceil(nb / bitsPerDigit) - 1;
  536.  
  537.     var b = biMultiplyByRadixPower(y, n - t);
  538.     while (biCompare(r, b) != -1) {
  539.         ++q.digits[n - t];
  540.         r = biSubtract(r, b);
  541.     }
  542.     for (var i = n; i > t; --i) {
  543.     var ri = (i >= r.digits.length) ? 0 : r.digits[i];
  544.     var ri1 = (i - 1 >= r.digits.length) ? 0 : r.digits[i - 1];
  545.     var ri2 = (i - 2 >= r.digits.length) ? 0 : r.digits[i - 2];
  546.     var yt = (t >= y.digits.length) ? 0 : y.digits[t];
  547.     var yt1 = (t - 1 >= y.digits.length) ? 0 : y.digits[t - 1];
  548.         if (ri == yt) {
  549.             q.digits[i - t - 1] = maxDigitVal;
  550.         } else {
  551.             q.digits[i - t - 1] = Math.floor((ri * biRadix + ri1) / yt);
  552.         }
  553.  
  554.         var c1 = q.digits[i - t - 1] * ((yt * biRadix) + yt1);
  555.         var c2 = (ri * biRadixSquared) + ((ri1 * biRadix) + ri2);
  556.         while (c1 > c2) {
  557.             --q.digits[i - t - 1];
  558.             c1 = q.digits[i - t - 1] * ((yt * biRadix) | yt1);
  559.             c2 = (ri * biRadix * biRadix) + ((ri1 * biRadix) + ri2);
  560.         }
  561.  
  562.         b = biMultiplyByRadixPower(y, i - t - 1);
  563.         r = biSubtract(r, biMultiplyDigit(b, q.digits[i - t - 1]));
  564.         if (r.isNeg) {
  565.             r = biAdd(r, b);
  566.             --q.digits[i - t - 1];
  567.         }
  568.     }
  569.     r = biShiftRight(r, lambda);
  570.     // Fiddle with the signs and stuff to make sure that 0 <= r < y.
  571.     q.isNeg = x.isNeg != origYIsNeg;
  572.     if (x.isNeg) {
  573.         if (origYIsNeg) {
  574.             q = biAdd(q, bigOne);
  575.         } else {
  576.             q = biSubtract(q, bigOne);
  577.         }
  578.         y = biShiftRight(y, lambda);
  579.         r = biSubtract(y, r);
  580.     }
  581.     // Check for the unbelievably stupid degenerate case of r == -0.
  582.     if (r.digits[0] == 0 && biHighIndex(r) == 0) r.isNeg = false;
  583.  
  584.     return new Array(q, r);
  585. }
  586.  
  587. function biDivide(x, y)
  588. {
  589.     return biDivideModulo(x, y)[0];
  590. }
  591.  
  592. function biModulo(x, y)
  593. {
  594.     return biDivideModulo(x, y)[1];
  595. }
  596.  
  597. function biMultiplyMod(x, y, m)
  598. {
  599.     return biModulo(biMultiply(x, y), m);
  600. }
  601.  
  602. function biPow(x, y)
  603. {
  604.     var result = bigOne;
  605.     var a = x;
  606.  
  607.     while (true) {
  608.         if ((y & 1) != 0) result = biMultiply(result, a);
  609.         y >>= 1;
  610.         if (y == 0) break;
  611.         a = biMultiply(a, a);
  612.     }
  613.     return result;
  614. }
  615.  
  616. function biPowMod(x, y, m)
  617. {
  618.     var result = bigOne;
  619.     var a = x;
  620.     var k = y;
  621.     while (true) {
  622.         if ((k.digits[0] & 1) != 0) result = biMultiplyMod(result, a, m);
  623.         k = biShiftRight(k, 1);
  624.         if (k.digits[0] == 0 && biHighIndex(k) == 0) break;
  625.         a = biMultiplyMod(a, a, m);
  626.     }
  627.     return result;
  628. }
  629.  
  630.